Periodic Vacuum Model and the QCD Instanton Vacuum 
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We formulate a matrix model for the 9 vacuum that embodies the n-point moments of the even and 
odd topological charge densities. The partition function is 2% periodic only after resumming over 
all winding modes. We apply the same analysis to the quenched instanton vacuum at finite 9 and 
derive saddle-point results for the vacuum energy density that are manifestly 2n periodic whatever 
9. We briefly discuss the consequences of unquenching the fermions on our results. 



PACS number(s); 12.38.Aw, ll.15.Tk, 11.30.-j. 

1. The non-abelian nature of QCD allows for the 
existence of classical, finite action, large gauge configura- 
tions that interpolate continuously between topologically 
inequivalent vacua. The canonically quantized vacuum 
state is quasi-periodic under large gauge transformations, 
with 9 playing the role of the quasi-momentum. States 
with different 9 correspond to different worlds. The dual 
to 9 is the topological charge and large gauge configura- 
tions of unit topological charge are instantons ||l ,^ . 

In the past years, the instanton approach to the QCD 
vacuum and its low-lying excitations has received consid- 
erable vindication both theoretically and numerically . 
In the instanton vacuum, 9 plays the role of an imaginary 
chemical potential, weighting incoherently the gauge con- 
figurations of net topological charge. As a result, most 
importance sampling arguments fail at finite 9. For small 
9, theoretical arguments based on the instanton vacuum 
yield a nontrivial 9 dependence of the vacuum partition 
function Q]. However, for large 9 the results lack explicit 
2-K periodicity, a point first addressed in 

In this study, we show how to construct a partition 
function for the quenched instanton vacuum that is man- 
ifestly 27r periodic whatever 9, following recent and gen- 
eral suggestions in Q. In section 2, we illustrate this 
approach by showing how the 27r periodicity is enforced 
in a matrix model of the QCD partition function, mak- 
ing the arguments of transparent. In section 3, we 
extend these arguments to the instanton vacuum model. 
In the quenched case, we derive explicit expressions for 
the instanton vacuum energy whatever 9. For small 9 our 
results reproduce those in [Q. In section 4, we extend our 
analysis to the unquenched case both in the matrix model 
and instantons. Our conclusions are in section 5. 



2. As finite 9 angle simulations for QCD are not fea- 
sible on the lattice, a study of the 6'-dependence is more 
accessible through a model such as the matrix model that 
can mock up the essentials of the 9 vacuum. The matrix 
model for the 9 vacuum was studied in H] with gaussian 
weights for the even and odd topological charge densities. 
Results from that model calculation directly address the 
levelling of the free energy as a function of 9, finite-size 
effects and effects of precision on the numerical analysis. 



However, we can inject more information into the model 
by using the low energy theorems to ensemble average 
with non-gaussian weights, as follows. In QCD, the con- 
nected 2-point correlation function for the scalar glueball 
obeys the low-energy theorem 
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:{T*^GG{x)^GGm = 



-4(|-GG(0)) , 



(1) 



where a is the strong coupling constant, and [3 is the 
1-loop beta function. The 2-point correlation function 
for the pseudoscalar and scalar-pseudoscalar mixing obey 
similar low-energy theorems 



d^x {T* — GG{x) 
Stt 



-GG(0)) 



e{^GGm 



(2) 



d'x{T*^GG{x)^GGm 



2e(£GG(0)), (3) 



where ^ q/p = A/b for the quenched case with 
h = -2iT(3/a^ = UNc/S. We note that in non-SUSY 
Yang-Mills theory, the precise values of p, q beyond one- 
loop are still debated |11|, therefore we refer to them 
generally as p and q in what follows. Higher correlation 
functions using the scalar and pseudoscalar sources can 
be defined using similar arguments |p^ . 

The matrix ensemble used to carry out the averaging 
T,„, which relate directly to the low energy 



has sizes n+ x 
theorems, thanks to the identification 

n± K [ d^x ( GG(x) ± iGGh 

lOTT J \ 



(4) 



which is a standard identification in matrix models of the 
instanton vacuum where n± counts the number of right- 
handed and left-handed zero modes (restricted to zero 
dimension) respectively. Then the low-energy theorems 
Eqs.(^,(||),(||) translate to conditions on the n- moments 



(5) 



Given the connected n- moments we can construct 
the associated generating function as 

W{J-, J+) = ^(n_)cxp(4J_) + ^(n+)exp(4J+) , (6) 
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where J-,J^ are pertinent sources. At finite vacuum 
angle 9, the generating function (^) changes. Since the 
low-energy theorems (0) also hold in the 9 vacuum, the 
9 dependence amounts to a change in the sources J± . 
This can be seen as follows. The 9 parameter enters as 
the coefficient of the topological charge 



Cg = frac9alphal6'K / d^x GG{x) 



(7) 



In terms of the matrix model, it is a coefficient of the 
topological susceptibility x ~ (n+— n_). When exponen- 
tiated to obtain the generating function, it thus serves to 
shift the sources for both ri+ and ri_, i.e. J_|_ and J_ 
as J± J± ±i£,9/2. Since ^ ^ ^/Nc, the dependence 
of the partition function on 9 appears as cos {9/Nc) (see 
below), which is clearly not periodic in 9. 



Z{9)^e^o^'^+'-^-^ =exp 



\ -(n+ + n)cos(-9)] , 
\1 P J 



(8) 



where we have used < — >= 0. This problem pre- 
cisely was addressed in A solution of the U{1) prob- 
lem in QCD requires that the 9 dependence does appear 
in this form However, physical quantities should be 
2TT periodic in 9 since it is a vacuum angle. This problem 
was identified by the authors of |^ by noting that the infi- 
nite volume limit taken in calculating physical quantities 
like the free energy does not commute with the 2tt shift 
in 9, and that it can be cured by imposing the constraint 
that the topological charge can take only integer values. 
This explicit global equation of constraint does not affect 
the physics extracted from the partition function defined 
in such a way. Therefore, to enforce manifest 27r period- 
icity in 9 in our model, we follow the arguments in 
and use the identity 



5(-(n+ — n_) — m) 
P 



m— — oo 



+ CXD 

E 

k— — oo 



ik — (n+ —71- ) 



(9) 



Note that this condition involves the integers p and q de- 
fined previously. All we are doing is introducing a diver- 
gent normalization ^ 1 to the definition of the partition 
function. However, this is not a trivial step since n+ and 
n_ are restricted to certain values satisfying the above 
condition. With this explicit summation, the partition 
fimction becomes 



kq. 



Z{9)^ exp ( -(n+-f n_)cos(— 6'-f 27r— ) ) , 

k— — oo ^ ^ 



(10) 



We may rewrite this as 



p-i 



1=0 



Z(6') = ^exp -(n++n_)cos( — 6i + 27r-)) , (11 



using the 2tt periodicity of the cosine. The vacuum en- 
ergy at finite 9 follows from E = W/V = —\nZ/V. 
Specifically, 



-1 /^"^ 

E{0) = ^lim —In i^^cxpiVEiO) cos{—6 
\i=o ^ 



(12) 



which is exactly the closed form version of the vacuum 
energy discussed in |^ using different methods. Here, 
E{0) is the vacuum energy at = 0. In fact, this is 
the functional form that must be obtained if the holo- 
morphic combinations GG ± i GG are selected to obtain 
the moments. In other words, the cosine is generic. The 
non-comniutativity of the thermodynamical limit and the 
shift 9 6 + 2tt is a, consequence of the multi-valuedness 
over I, as pointed out first in This means that peri- 
odicity in 9 can be reconciled with the thermodynamical 
limit only if we select the right values of In simple 
terms, we encounter different "Reimann sheets" in mak- 
ing the 2tt shift so that the entire branch structure must 
be retained to account for the periodicity in 9. Armed 
with a partition function for the matrix model that has 
the correct periodicity in 0, we can comment on the fea- 
tures that would arise in a numerical simulation of such 
a model. We stress that a matrix model is not intended 
to represent the dynamics of the QCD vacuum. Rather, 
we have written down a partition function within a toy 
model that is amenable to numerical simulations which 
would yield the qualitative dependence of the free energy 
on 9. 

In anticipation of possible numerical simulations of the 
partition function at finite 9, it is useful to analyze the 
partial partition functions for configurations with N = 
maxn-t. Specifically, 
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Zni(9) = cMVE(0) cos(--9 + )) 

P P 

i,p)r{N + i,i^) 



r{N - 



Nl 



Nl 



where 



and 



II = (n+)- exp(-t-9 H ) 

q P P 



, \P /-In 27:1. 

(n„)-exp(z-6l-)- ) 

q p p 



(13) 



(14) 



(15) 



Here, T is the incomplete gamma function. Note that for 
— > cx) we recover the partition function (|ll|) . Recently, 
it was argued by a number of authors [g 12 1^ that a 
lattice assessment of the vacuum energy at finite 9 in 
CP" models and Yang-Mills theory show an unexpected 
levelling at large 9 with consequences on the strong CP 
problem. However, this observation may be a numeri- 
cal artifact |§,M. In our case, this issue can be readily 
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investigated by noting that ( p^ approximates the full 
partition function at finite but large N. Numerically, 
corrections in start to be appreciable when 



In (q/pn^) 



71*- cos{q9/p) < 1 



(16) 



where — N/V is the maximum density of winding 
modes. For simplicity, we have specialized to the I — 
branch for which the ground state occurs for < 6* < 
For fixed p and q, the deviations between the results for 
N = oo and N finite decrease with increasing n*. For 
fixed q and n*, the deviations decrease with increasing 
p (recall that p ~ Nc). The changes in the p,q,n^ pa- 
rameters in the finite sum (13) can significantly alter the 
9 dependence of the vacuum energy (|jl^. These ob- 
servations are useful for the lattice calculations carried 
in ||l^,|l^. We now move to applying this periodicity 
procedure for 9 to the interacting instanton ensemble. 



3. The prescription for ensuring 2tt periodicity for the 
model vacuum partition function discussed above, can be 
readily extended to a weakly interacting ensemble of in- 
stantons and anti-instantons For that, we consider 

the vacuum partition function at finite 9 (vacuum angle) 
and /i (chemical potential for winding modes) from H] 



Z{fi,9) 



N± 



(17) 



The partial sums are given by 



Zn± — 



N+l 



n 



(18) 



where (± are the pertinent fugacities yj. The instantons 
are characterized by a set of collective coordinates (i?), 
such as their position, color orientation and size (p). For 
simplicity, we choose the relative two-body action in the 
schematic form 

S(i?l,i?2) ~ -^712 iPlP2f , 

as suggested by variational calculations. Here 712 — js 
for like particles, and 712 = 7a for unlike ones. A is the 
inverse charge at the typical relative separation and V is 
the 4-volume. The 9 dependence of the vacuum energy 
is easy to write down for special values of r = ^a/lsi 
pointed out in [Q. For r = 1, (which is the sum ansatz 

of i) 



f{9) = [cos(0)]' 



(19) 



where f{9) expresses the 0-dependence of the partition 
function written in a generic form as 



Z{^,, 9) = exp(^(7V)o/(0)exp(^)) (20) 

and (A^)o is the mean instantonf/) and anti-instanton(/) 
number at zero vacuum angle 0. For r = 0, 



46* 

/(0)=cos(-), 



and finally, for r 
f{0) 



-1, 



261 

6-2 



(21) 



(22) 



The key point we notice is that for r = and r —1, 
the free energy is not 27r periodic in 9. The problem and 
it's solution are clear once we see that the free energy 
is extracted following a saddle-point analysis discussed 
in 1^]. It involves the relative ratio 



ln(iV+/7V_) 
2i 



which satisfies the saddle-point equation 



-arctan- 



(1 - r^)sin$ 



cos $ + rV 1 — r^sin^ $ 



(23) 



(24) 



A shift of 27r in 9 can be compensated by a similar shift 
in (f>, so 2tt periodicity seems implicit here. However, 
the free energy involves taking a thermodynamic limit 
which does not commute with the 2tt shift as explained 
before. This is why the free energy expressions are not 
27r periodic in 9 (except for r = 1 which is a special case 
discussed later). If the free energy, which is a physical 
quantity, has to be 27r periodic in 9, the explicit but un- 
conventional sum over winding modes has to be enforced 
as before, and the modified partition function is 

+00 

k— — oo N± 

(25) 

The fixed N± partition function can be varied to deter- 
mine the best size distribution. Evaluation of the parti- 
tion function requires an expression for the average size 
as well, and these can be determined from consistency 
equations. The full partition function can then be cal- 
culated in a saddle-point analysis. We refer the reader 
to Q for a further and explicit discussion and state here 
that the new saddle point equation becomes 



<3> H ■ — arctan 



(1 - r^)sin<i> 



cos $ + r\/ 1 - r^sin^ $ 



= 9 + 2TTk 



(26) 



The infinite sum over k looks superfluous but it is not. 
It implies an infinite number of saddle points, but as 
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before, this may be rewritten as a principal solution in 
a branch /, plus an infinite number of copies, that can 
be dropped because they all yield the same minimum. 
This also explains why the summation over I is finite ( 
to p — 1 ). The value of I to choose depends on the 
range of 9 and the integers q and p. These observations 
are novel and absent in [Q. If we ask the question of 
2TT periodicity of 9 after performing the thermodynamic 
limit, it is essential that this limit is taken properly, that 
is, with all branches described by I taken into account. 
When the thermodynamic limit is taken properly, the 
character of the vacuum energy at finite 9 is changed 
and is explicitly periodic as we now show. For r — 0, 



f{9) = cos — + ) , 

p 



(27) 



where I may be any integer between and p — 1, depend- 
ing on which range 9 lies in. To see the 27r periodicity, one 
has to retain all appropriate branches of 1 corresponding 
to the different ranges of 9 as it moves through 2tt. Then 
we observe that the soutions are merely relabelled but 
the entire set in / remains the same. This is the subtle 
way in which periodicity in 9 is enforced. Finally, for 



— 1, we get 



f{0) = 



cos 



29 

b-2 



2Trl 
P 



2(b-2) 



(28) 



Here too, I lies between and p — I, but the ranges for 
these choices is different because of the different cosine 
dependence. This is a limiting case since a stable ensem- 
ble does not exist for r < — 1 . Eqns (p^) , (p7|) , (|2^) are new 



and flow from the modified partition function Eq.(25). 
The case r = is special since in this case, the expression 
is periodic without the 27rfc which is indeed superfluous. 
This seems to indicate that the free energy depends only 
on the total number of pseudoparticles in the system, and 
is insensitive to the difference in the number of instan- 
tons and anti-instantons. This would be natural since the 
sum ansatz sets the II interaction equal to the // inter- 
action. In conclusion, we see that the results of |Q admit 
a generalization in our approach which makes them man- 
ifestly periodic in the same way as in the model discussed 
above. Lattice cooling measurements with rich instanton 
content should be useful in assessing the 9 dependence of 
the vacuum energy along the present lines, including the 
value of r. 



4. In so far, both the model and instanton analy- 
sis of the energy of the 9 vacuum was carried without 
fermions (quenched). The addition of fermions modify 
substantially our results. Given the similarities between 
the model and the instanton analysis, we will show in 
this section how our model calculation is affected by the 
introduction of fermions. This is of course a simplified 



version of a full fermion analysis in an instanton ensem- 
ble, but should be useful for an unquenched lattice sim- 
ulation nevertheless. We will derive analytic results for 
the vacuum energy as a function of 9 for small fermion 
masses in the cases Nf = 1,2. From this, the topologi- 
cal susceptibility can be determined and our model can 
be confronted with experimental cooling simulations on 
a lattice. The addition of fermions {Nf flavors) modifies 
the model partition function as follows |p|,p^ 



Z{9,Nf) 



(ndct( 



w 



)) (29) 



irrije 



where W is a complex hermitian asymmetric matrix of 
dimension fc+ x fc_ (not the same as n+, n_). The reason 
that k± is not the same as n± has to do with the chiral 
anomaly, as we will see. For simplicity, the measure on 
W is dictated by the maximum entropy principle, hence 
gaussian. The quark mass matrix is introduced in a form 
consistent with the chiral anomaly 



9 — > 6* arg det M 



(30) 



To perform the disorder average successfully, we must 
make a modification in the process of bosonization, that 
will determine the weight to be used. In the presence 
of fermions, the explicit 27r periodicity is enforced by re- 
quiring that the matrix asymmetry satisfies 



fc+ — fc_ = m 



(31) 



with m fixed by the additional identity (|9|). This con- 
straint implies that the number of zero modes is not the 
same as the number of winding modes for the unquenched 
case, and is essential to preserve the form of the chiral 
anomaly. Now, bosonization may be performed in a man- 
ner consistent with the anomaly, and the unquenched 
partition function reads 

1, ^ 



p-i 
x^e 



Vi5(O)cos(^(e+0)+27ri) 



(32) 



where k = fc-|. -I- fc_ and cj) = argdetM. Once again, the 
periodicity in 9 is realized over the set of all branches of 
the function, so we pick a particular value of I for each 
range of 9 + (p. For example, when < {9 + (f)) < Tr/q, 
the I = branch dominates in the thermodynamic limit. 
A saddle-point approximation to the partition function 
may be carried out with P = diag pje~^^'^'~^^'^'\ (where 
i labels flavor) to yield the unquenched vacuum energy. 

Considering the range of 9 for which I — dominates 
and for small masses rrij <g; pj , we obtain for the un- 
quenched vacuum energy 
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Nf 



i=i 



m = T E - P] + 2^ cos 0,] 



+E(Q) cos(~ -e + ^(f>.j + !^sm(f>A] (33) 
P P P] J 

to 0{m^), where we have taken k to be peaked around 
fc*. Thus, we have written down a distribution in fc+ + fc_ 
apart from the natural one in fc+ — fc_ . The saddle points 
in the p's decouple to give 



Pj = 1- {^) cos e + 0{m^) 



(34) 



and to the same order, the saddle points in the cj)'s is 



= Y,<j,,+O{m) . 



(35) 



It is easy to show that 



TTii sin 01 = TO2 sin 02 = ...tojvj. sin (jj^^ . (36) 

These approximate expressions were obtained in the large 
Nc limit as well Note that the topological suscepti- 
bility X for equal quark masses follows from 

£^,„{9) = -E{0) - K mNf cos{e/Nf) + 0{m^) (37) 

in the form 

a well-known result . 

In general, for the 1th branch 

(2;-l)^<(0-^0,)<(2/ + l)^ , (39) 

and unequal masses, we have 



E{e) = -£;(0)cos( -^d+^Y 
\ p p ^ — ' 



N, 



27r- 



(40) 



For Nf = 1, we have 



f ^* (If -l]sme + 0{m^) , (41) 



for which the subtracted vacuum energy reads 

A£{e) = k^m{l-cos9) . (42) 

Although the branch-structure affects the position of the 
saddle point, it does not produce a physical change in the 



vacuum energy. At small m, for any 9, there is no cusp. 
There would be a cusp at = tt if to were large but 
the saddle point analysis is valid only for small m. For 
Nf — 2, we have 



sin 01^2 = ± 



TO2,i sm( 



\/'ml + : 



2toiTO2 cos ( 



(43) 



with the subtracted vacuum energy 



A£{9) 



I 7711 + 7772 I - \Jml + ml + 



27771 TO2 cos 9 ] . 

(44) 
TT, which 



For TOi = 7772, the free energy has a cusp at 9 
indicates a spontaneous breaking of strong CP. 

Most of the above results apply to an ensemble of 
weakly interacting instantons and anti-instantons as dis- 
cussed above, since the pair-interaction is "screened" by 
the quarks through the fermion determinant. The mul- 
tivalued structure of the unquenched partition function 
should be similar to what is shown in this section. 



5. In conclusion, following on recent suggestions 
we have shown how a multivalued approach to the 9 
vacuum can be displayed in a simple model, deriving 
closed form results for the vacuum energy that are man- 
ifestly periodic whatever 9. This model is particularly 
useful for testing the pertinence of numerical simula- 
tions in a field theory (CP" and Yang-Mills) at finite 9. 
We have used the insights generated by the analytically 
tractable model to derive the vacuum energy at finite 
of a quenched and dilute instanton anti-instanton en- 
semble, that is manifestly 27r periodic, generalizing ear- 
lier results 0. While direct lattice calculations at finite 
9 are not feasible due to the complex integration mea- 
sure, quenched and cooled simulations can test the en- 
ergy dependence on 9 via a measurement of the topolog- 
ical charge. We have shown how the fermions affect the 
model calculation, and suggest that the same generic re- 
sults apply to the unquenched instanton calculation. Our 
analysis can be improved, by considering £, = q/p beyond 
the one-loop result, still a point of discussion pi]]. 
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